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Abstract 



cs . .... 

J> , We derive the time-evolution equation that describes the Brownian motion of labeled individual 

0^ ' tracer particles in a simple model atomic liquid (i.e., a system of N particles whose motion is gov- 

00 ' 

erned by Newton's second law, and interacting through spherically symmetric pairwise potentials). 

en . 

^^ . We base our derivation on the generalized Langevin equation formalism, and find that the resulting 

time evolution equation is formally identical to the generalized Langevin equation that describes 
the Brownian motion of individual tracer particles in a colloidal suspension in the absence of hydro- 
dynamic interactions. This formal dynamic equivalence implies the long-time indistinguishability 

a ■ 

of some dynamic properties of both systems, such as their mean squared displacement, upon a 
well-defined time scaling. This prediction is tested here by comparing the results of molecular and 
Brownian dynamics simulations performed on the hard sphere system. 



I. INTRODUCTION. 

It is well known that under some circumstances the phenomenology of atomic liquids finds 
an almost perfect correspondence in the phenomenology of colloidal fluids \MiM- This seems 
to be particularly true regarding the rather complex dynamic behavior of these systems as 



01?]. 



they approach the glass transition [5H7J. Although it is clear that this analogy has some 
fundamental limitations (such as the presence of many-body hydrodynamic interactions in 
colloidal systems), one can be confident, for example, that the phase behavior of colloidal 
and atomic systems with identical interaction potentials will, of course, be the same. Thus, 
if one approaches this problem with a dynamic simulation technique, one is confident that 
the equilibrium phase diagram of a specific system (say a Lennard- Jones liquid) will be 
independent of the simulation technique employed in its determination (either molecular or 
Brownian dynamics) [8J]. Time-dependent and dynamic properties, on the other hand, are 
expected in general to depend on the specific microscopic transport mechanisms. Never- 
theless, some features associated with the collective, long-time behavior of the system also 
seem to be rather insensitive to the microscopic short-time dynamics. For example, it has 
been suspected, and partially corroborated, that for a given model system (i.e., same pair 
potential) standard molecular dynamics will lead to essentially the same dynamic arrest 
scenario as Brownian dynamics [5H3]. Determining the range of validity of this dynamic 
analogy, however, continues to be a relevant topic in the study of the dynamics of liquids. 

From the theoretical side, for example, one would like to have a unified description of 
the macroscopic dynamics of both, colloidal and atomic liquids, which explicitly exhibits 
the origin of the similarities and differences in their macroscopic dynamics. This topic has 
jeen addressed in the framework of the mode coupling theory of the ideal glass transition 
9) , originally developed for Newtonian liquids, but also adapted to Brownian systems. Such 
attention was focussed on the similarity of the long-time dynamics of Newtonian and Brow- 



nian systems in the neighborhood of the glass transition 1QJ. A number of issues, however, 
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remain open |5j], one important question referring to the validity and limitations of this 
long-time similarity under general conditions, such as those involving ordinary thermody- 
namically stable fluids, and not necessarily associated with the glass transition. In this sense, 
one possible general framework for such theoretical analysis is the concept of the generalized 
Langevin equation (GLE) [111. Il2|. 



The GLE formalism describes the dynamics of the thermal fluctuations 8di(t) (= 
di(t) — a1 q ) of the instantaneous value of the macroscopic variables a,i(t) (i = 1,2, ...,z/), 
around its equilibrium value a^ 9 , and has the structure of the most general linear stochas- 
tic equation with additive noise for the vector 8&(t) = [8a\(t), Sa^t), ..., 5a u (t)y (with the 
dagger indicating transpose). The GLE equation has been widely used in the description 
of thermal fluctuation phenomena in simple liquid systems, and Boon and Yip's textbook 
131 ] contains a detailed account of its early use to describe the dynamics of simple liq- 
uids. Although this stochastic equation is conventionally associated with the Mori-Zwanzig 



projection operator formalism 14J, |l5|, in reality its structure is not a consequence of the 



hamiltonian basis of Mori-Zwanzig's derivation; instead, it is essentially equivalent to the 



mathematical condition of stationarity 12]. 

Understood in the latter manner, the GLE formalism was first employed in Ref. [11 1 

to derive the equation of motion of an individual tracer particle in a colloidal suspension 

without hydrodynamic interactions. Such an equation reads 
r dv(t) 



Ni- 



di 



-C (s) v(£) + fW(£)- I dt'At(t-t')v{t')+F{t), 
Jo 



(1.1) 



where M is the mass and v(t) the velocity of the tracer particle, while Q s > is the friction 
coefficient caused by the frictional resistance of the supporting solvent and f (*' (t) the associ- 
ated random force. The memory term involving the time-dependent friction function A((t), 
and its associated random force F(t), are the friction and fluctuating forces that originate in 
the time-evolution of the cage of surrounding colloidal particles. Under well defined approx- 
imations, the exact result for the time-dependent friction function A£(£) derived in Ref. 11 ] 
was shown there to reduce to the following approximate expression in terms of the collective 
and self intermediate scattering functions (ISFs) F(k,t) and F$(k,t), 
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dk 
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F(k,t)F s {k,t). 
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In this equation T is the temperature, n the number concentration, and S(k) the static 
structure factor of the bulk suspension. This result, together with similarly general expres- 
sions for F(k, t) and F$(k, t) also derived within the GLE formalism 16], was later employed 



in the construction o 
cof 



oid dynamics [17 



the self-consistent generalized Langevin equation (SCGLE) theory of 



181 ] , eventually applied to the description of dynamic arrest phenomena 



19H21J. and more recently, to the construction o 



and aging of colloidal glass-forming liquids 
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23]. 



a first-principles theory of equilibration 



With the aim of investigating the relationship between the dynamics of atomic and Brow- 
nian liquids, in this work we apply the GLE formalism to derive the generalized Langevin 
equation that describes the motion of individual tracer particles in simple atomic liquids, 



thus extending to these systems the results of Ref. llj reviewed above. The most remark- 
able prediction of the derivation presented here is that the resulting stochastic equation for 
the velocity v(t) of the atomic tracer turns out to be formally identical to the colloidal case 
described by the two equations above, with the solvent friction coefficient C replaced by a 
kinetic (or 'Doppler') friction coefficient (° determined by kinetic-theoretical arguments. 

Since the concept of kinetic friction may involve a rather subtle use of otherwise simple and 
well established concepts, we start this paper in section II by providing a simple and intuitive 
description of the short-time random motion performed by an individual tracer particle in 
an atomic liquid, as a consequence of molecular collisions. Such description exhibits the fact 
that the resulting random motion must be described by the same stochastic mathematical 
model that describes the Brownian motion of a tracer particle in a colloidal fluid. This is 
just the mathematical model underlying the ordinary Langevin equation (Eq. (II. ip above 
without the time- dependent friction term and its associated random force). Thus, also in the 
atomic case, a relaxation time To of the velocity, due to the friction force — C°v(£), defines the 
crossover from ballistic to diffusive motion. The fundamental difference lies in the physical 
origin of the friction force — C° v (^) and in the determination of the friction coefficient £°- In 
a Brownian liquid the friction force — C v(t) is caused by the supporting solvent; hence, C 
assumes its Stokes value. In contrast, as discussed in section II, in a Newtonian liquid the 
friction force — C° v (£) is not caused by any external material agent but by the unimpeded 
tendency to establish or restore, through molecular collisions, the equipartition of the energy 
available for distribution among the kinetic energy degrees of freedom of the system. Thus, 
the corresponding value of (° that emerges from these considerations is provided by Einstein's 
relation, with a kinetic-theoretically determined diffusion coefficient. 

The discussion of the Langevin equation for atomic liquids is continued in sections IIHI 
and IIVI where a more formal and complete derivation is provided. Thus, in section IIHI we 
show that the force on the tracer particle can be written as an integral of the divergence 
of the stress tensor n (r, t), whose kinetic component Hk (r, t) is the origin of the Doppler 
friction force —(°v(t) and its associated random force, whereas the term involving the con- 
figurational component lit/ ( r 5 1) describes the effects of the ordinary conservative direct 



forces (electrostatic, van der Waals, etc,) exerted by the surrounding particles. It is also 
seen that the latter effects enter additively in the Langevin equation of the atomic liquid, 
as an integral term that is linear in the instantaneous local density of the surrounding par- 
ticles. In section IIVI we demonstrate that this linear coupling of the motion of the tracer 
with the local density of the surrounding particles leads to a time-dependent configurational 
friction term and to its corresponding random force, thus resulting in a generalized langevin 
equation for the velocity of a tracer particle, which turns out to be formally identical to its 
colloidal counterpart in Eq. (II. ip . 

The main predictions of the resulting generalized Langevin equation for atomic liquids are 
then discussed in Section |V] These include well defined scaling rules that exhibit the identity 
between the long-time dynamics of atomic and colloidal liquids. There we test these scalings 
by comparing the simulation results for a given model system (the hard sphere fluid) using 
both, molecular dynamics and Brownian dynamics simulations. The last section summarizes 
the most relevant conclusions, and discusses some limitations and potential applications of 
the results of the paper. 

II. BALLISTIC AND DIFFUSIVE REGIMES IN A SIMPLE ATOMIC LIQUID. 

Let us consider a simple atomic fluid, formed by N spherical particles in a volume V 
whose microscopic dynamics is described by Newton's equations, 

M^ = X>,-(t), (* = 1,2,...,A0, (2.1) 

i¥=i 

where M is the mass and Vj(t) = dri(t)/dt the velocity of the ith particle at position Ti(t), 
and in which the interactions between the particles are represented by a sum of pairwise 
forces, with F^- = — Wit(|r$ — r^l) being the force exerted on particle % by particle j. Our 
general aim is to establish a connection between the microscopic dynamics described by 
these (Newton's) equations, and the equation that describes the random motion of any 
representative individual tracer particle of the liquid. 

In this section we discuss a simple and intuitive (albeit possibly subtle) physical picture, 
which plays a central role in our theoretical effort to establish such connection. We start by 
recalling the physical meaning of two fundamental concepts in the dynamics of an atomic 
fluid, namely, the mean free time, tq, and the mean free path, l , which represent the 



characteristic time- and length-scales at which the crossover from the short-time bal 



istic 



motion of the atoms to their long-time diffusive transport occur. It is well known [13|, |24| 
that for correlation times t much shorter than r , and for distances much shorter than l , all 
the particles move ballistically, so that for t <C r , the mean squared displacement (MSD) 
is given by < (Ar(£)) 2 > ~ 3t> 2 £ 2 , with v = (ksT/M)^ being the thermal velocity. 

For times t much longer than ro, each particle has undergone many collisions, and its 
motion can be represented as a sequence of (ballistic) random flights of mean length l Q and 
mean flight-time tq, traveled at a random velocity that has zero mean and covariance v 2 . 



From the theory of random flights, however, it is well-known 25| that the motion represented 
by such a sequence of random displacements, will become diffusive in the long-time limit. 
This means that it will be characterized by a mean squared displacement that, for t ^> To, 
will increase linearly with time, < (Ar(£)) 2 > « 6D°t. Furthermore, the corresponding 
diffusion coefficient D° will be given by D° = (7o) 2 / r o- 

Thus, we conclude that the MSD of a representative tracer particle will exhibit two well- 
defined limiting behaviors in two opposite time regimes, namely, it will be ballistic at short 
times, < (Ar(£)) 2 > ps 3w 2 t 2 for t <C r , and diffusive at long times, < (Ar(£)) 2 > rj 6D°t for 
t ^> Tq. The simplest mathematical model that provides a full description of the crossover of 
< (Ar(t)) 2 > from the first to the second of these two exact limits, is provided by a Gaussian 
stationary stochastic process, described by a linear stochastic equation with additive noise 



for the instantaneous velocity v(£) of the tracer particle, i.e., by [26| 

M^- = -C°v(t)+f°(t), (2.2) 

with f°(t) being a "purely random" (or "white") noise, i.e., a stationary and Gaussian 



stochastic process with zero mean (f°(t) = 0), uncorrelated with the initial value vo of the 
velocity fluctuations, and delta-correlated with itself. In fact, the stationarity condition is 
in reality equivalent to the fluctuation-dissipation relation between the random and the dis- 



sipative terms in Eq. (ET2"]) . f°(t)f°(t') = (v v )MC 2(5(t-t / ) ) where (v v ) is the stationary 
covariance of the initial velocities. Within the additional physical assumption that identifies 
the stationary state described by Eq. (12. 2ft with the thermodynamic equilibrium state, we 
have that this covariance is determined by the equipartition theorem, (voVo) = (A;_bT ' /M) I 
(with I being the 3x3 cartesian unit tensor). 



Clearly, this equation is formally identical to the ordinary Langevin equation 27| for the 



instantaneous velocity v(t) of a colloidal particle in a solvent, in which case, Stoke's result 
(^ = Sitr/a (with 77 being the viscosity of the solvent and a the diameter of the colloidal 
particle) provides an independent determination of the friction coefficient Q 8 ^ [28]. In the 
present case, however, there is no supporting solvent to produce friction, and hence, identify- 
ing the origin and determining the value of (° requires slightly more subtle arguments. The 
simplest manner to describe its physical origin may be found in Uhlenbeck and Ornstein's 



brief reference to the so-called Doppler friction 29J. These authors point out that any tracer 
particle colliding with the particles of a gas, whose size a is smaller than the mean free path 
Iq, will be subjected to Doppler friction, caused by the fact that "when the tracer particle is 
moving, say to the right, will be hit by more molecules from the right than from the left" . In 
the following section we shall provide a more formal derivation of this kinetic friction effect. 
At this point, however, we provide simple arguments for its quantitative determination. 

To determine (° in the present case we first notice that the mathematical solution of Eq. 
(f2T2l for the MSD is such that at long times, < (Ar(£)) 2 > « 6(k B T/(°)t. This determines 
a self-diffusion coefficient D° in terms of (° through Einstein's relation, D° = fcsT/C - Thus, 
we must only determine either (° or D°. In our case, we write Einstein's relation as 

C° = k B T/D°, (2.3) 

and determine D° independently, borrowing the arguments developed in the elementary 



kinetic theory of gases [28|. For this, we recall that D° = (Zo) 2 / r O) which, since Zo/ r o = v o, 
can be written as D° = VqIq. We then estimate the mean free path l to be given by 
lo ~ l/ncr 2 , with n = N/V and with a being the collision diameter of the particles. Thus, 
we must have that D° ~ yksT ' jM /(na 2 ). In fact, the rigorous value of D° is [30|] 



3 fk B T\ 1/2 1 
7^ V~^T/ ™ 2 ' 



_ o li B 



^^ir —,■ (2-4) 



This expression, together with Einstein's relation above, determines the value of the kinetic 
friction coefficient (° of an atomic fluid. 

Let us now discuss some of the implications of the atomic Langevin equation defined by 
Eqs. ( 12. 2p . ( 12. 3p . and ( 12. 4p . Since the Langevin equation itself is mathematically identical 
to the ordinary (i.e., colloidal) Langevin equation, its solution is also formally the same. For 
example, from Eq. (12.2J) . and the assumed properties of f°(t) one can evaluate the velocity 



auto-correlation function (VAF) 

V(t) =< v(t) ■ v(0) > /3, (2.5) 



(2.6) 



with the result 


V(t) = v 2 e~ t ^, 


with 






T S = M/C 


being the velocity relaxation 


time. 


The MSD, normalized as 






W(t) =< (Ar(t)) 2 > /6, 


is related with the VAF by means of the exact relationship 




W(t)= f (t-t')V(t')dt' 
Jo 



(2.7) 



(2.8) 



(2.9) 

or, in terms of the Laplace transforms (LT) W(z) and V(z), as 

W(z) = V(z)/z 2 . (2.10) 

This exact relationship can be written, for the particular form of the VAF in Eq. (12.61) . as 
the following differential equation, 

r s d ^ + Wif) = D\ (2.11) 



whose solution reads 



W(t) = D°r s 



t 


t 

- 1 + e T s 


TS 





(2.12) 



This expression interpolates W(t) between its corresponding short- and long-time asymptotic 
limits, 

W(t) K-vlt 2 , fort<r (2.13) 

and 

W(t) w D°t, for t > T . (2.14) 

In addition, it also exhibits the fact that the crossover from ballistic to diffusive motion is 
most naturally described using t s as the unit of time. 



Let us notice, however, that for atomic liquids the relaxation time Ts is identical to the 
mean free time r , since t s = M/(° = (k B T/(°)(M/k B T) = D°/v$ = (^/r )(r // ) 2 = r . 
Thus, for an atomic liquid the mean free time is the most natural time unit, and the mean 
free path l = vqt the most natural length unit since, for an atomic liquid, Eq. ( 12.1 2p can 
be rewritten in terms of the scaled time t* = £/to and the scaled MSD w(t*) = W(t)/lQ as 



w (t*) =[t*-l + e - **] . 



(2.15) 



To illustrate the validity of this result, in Fig. 1 we compare it with the molecular dynamics 
simulation data (solid circles) for w(t*) in a fluid of hard spheres of diameter a at a small 
but finite volume fraction <p = Tcna 3 /G, namely, at <fi = 0.1. Clearly, the scaled solution 
(I2.15P of the atomic Langevin equation (solid curve) lies very close to the simulation data 
of the MSD. 
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FIG. 1: Scaled mean squared displacement [W(t)/Zjj] as a function of the scaled time [t/rg], with Iq 
and ts = tq being the mean free path and mean free time, of a hard-sphere fluid at volume fractions 
cf) = 0.1 (solid circles) and <p = 0.5 (empty circles) generated by molecular dynamics simulations of 
soft sphere systems using the velocity Ver let's algorithm [8] and the soft-to-hard-sphere dynamic 
correspondence of Ref. [3j|. The solid line corresponds to the low density limit, Eq. (]2.15p . and 
the dashed line corresponds to the expression in Eq. (|2.17p with D* = 0.099. 



The dimensionless quantities w(t*) and t* serve to highlight an important scaling that 
derives directly from the atomic Langevin equation in Eq. (12. 2p . At the same time, however, 
these scaled variables hide the specificity of the actual magnitudes involved in concrete real 



physical situations. For example, let us think of a typical atomic liquid, such as argon 
(M rj 40 amu, a ~ 3. 8 A), and rewrite the expression for D° in Eq. (j2.4p as 



D , ^ 



160 



(T 



M 



(2.16) 



Then, for T rj 300°^ and <p ~ 0.1, representative of a moderately dilute gas at room 
temperature, we find that D° rj 1.0 x 10 _7 m 2 /s, Zq ~ 4.2A, and 75 = To rj 1.7ps. Let us now 
compare this crossover timescale Ts with the structural relaxation timescale 77 = d 2 /D°, the 
time it would take any particle to diffuse one interparticle mean distance d = n^ 1 ^ with the 
diffusion coefficient D°. For the conditions under consideration we find that 77 rj 4.16ps, so 
that tj is only about twice longer than t$. This means that there is nothing like a time-scale 
separation between the crossover from ballistic to diffusive motion and the crossover from 
free to collective relaxation (i.e., from the so-called /3-processes to the slower a-processes). 
Although this timescale degeneration disappears gradually as the density increases, giving 
rise to the full separation of time scales characteristic of the approach to the glass transition, 
this degeneration of time scales is a characteristic feature of atomic dynamics throughout 
the whole stable liquid phase of these simple systems. To illustrate this, notice that taking 
T rj 83° K and rj 0.5, representative of the freezing conditions of Argon, leads to D° rj 
0.11 x 10~ 7 m 2 /s, l rj 0.84A, ts = t rj 0.64ps, and 77 RJ 13.5ps., so that now 77 rj 20ts. 

This situation must now be compared with its corresponding colloidal analog. In this case, 
the velocity relaxation time t$ = M/(° is no longer identical to the mean free time. Instead, 
it is determined by the mass M of the particle and by the solvent friction coefficient Q"\ 
now given by its Stokes value Q") = 3nr]a, with r\ being the viscosity of the solvent. Thus, 
consider a micron-sized colloidal particle in water at room temperature and with the same 
mass density as water itself, so that o rj l.O/i, M rj 0.52 x 10~ 15 kg, and (^ rj 0.94 x 10 -8 
kg/s (taking rj rj 10~ 3 kg/m-s for the viscosity of water). This leads to the estimate D° rj 
4.4 x 10 _13 m 2 /s for the (^-independent) Einstein diffusion coefficient D° = kBT/Q s ' . The 
corresponding estimate of r$ = M/^ is, then, t$ ~ 5.5 x 10 _8 s, whereas the value of 
77 = d 2 /D° (for = 0.1) is 77 ~ 6.87s. We thus immediately notice a really dramatic 
separation of timescales, expressed by the fact that now 77 ~ 10 8 7"s. As a consequence, 
in a colloidal suspension, the crossover from ballistic to diffusive motion (occurring at very 
early times t rj t$) and the crossover from free to collective difusion (occurring at times 
t ~ 77), are separated by nearly 8 decades. Thus, if one were interested in observing the 

10 



ballistic-to-diffusive crossover in our illustrative colloidal suspension, one would have to 
measure W(t) in the time window corresponding to Fig. 1, i.e, for times centered in the 
regime t « t$ ~ 10~ 8 s. This crossover, however, will be completely shifted to unobservable, 
extremely short times, when probing only the structural relaxation of the colloidal liquid for 
times t in the milisecond range and above (i.e., t rj tj). For this reason in the description of 
the dynamics of suspensions the inertial term of the Langevin equation is normally neglected, 
M-^p ~ 0, leading to the overdamped Langevin equation, —(^dt + f^W = 0> whose 
solution for the MSD is W(t) ~ D°t. This, however, is just the diffusive limit in Eq. (j2.14p . 

The main conclusion of these illustrative estimates is, thus, that in colloidal liquids this 
wide separation of timescales is always present, whereas it is virtually nonexistent in an 
atomic liquid. As mentioned above, however, this disparity between atomic and colloidal 
liquids will disappear in the metastable regime due to the dominant effect of interparticle 
interactions in the determination of the structural relaxation processes. To introduce the 
discussion of these effects, in Fig. 1 we have also included the molecular dynamics simulation 
results for the MSD of the hard-sphere liquid at freezing conditions, <ft — 0.5 (empty circles). 
These results illustrate the deviations from the ideal behavior described by Eq. ( 12.151) . 
which originate from the structural effects of the interparticle interactions (very moderate, 
and almost imperceptible in the results for = 0.1). The comparison between the results 
for the MSD at these two volume fractions clearly show that in both cases the ballistic 
short-time behavior (i.e., for times t < r ) is accurately described by Eq. (j2.15p . Beyond 
the crossover time To, however, the results for <fi = 0.1 exhibit only the kinetic effects of the 
interparticle interactions, whereas the results for = 0.5 exhibit the combination of both, 
the kinetic and the structural effects of the interparticle forces. 

One of the results of the following sections will be a simple expression for W(t) that 
preserves the short-time limit W(t) ~ \v 2 t 2 , but replaces the ideal long-time limit W(t) ~ 
D°t by the correct one, W(t) ~ -D^t, where Dl is the long-time self-diffusion coefficient (see 
end of Sec. HVl) . This expression reads 

w (t*) = D* 2 [t*/D* - 1 + e -'* /D *] , (2.17) 

with the parameter D* = Dl/D° being the ratio of the long-time to the short-time self- 
diffusion coefficients. In fact, it is the solution of 

*p. + ap. = e, (2.i8) 

dt* D* K ' 

11 



which extends Eq. (12. lip to finite densities. The dashed line in Fig. 1 corresponds to this 
expression with the value D*=0.099. The comparison with the simulation data indicates 
that Eq. (I2.17P provides an accurate representation of W(t) at short and at long times 
compared with r . In this comparison, the parameter D* was employed as a mere fitting 
parameter. The idea, however, is to build a first-principles theory that predicts the value of 
D* and the behavior of W(t) in the intermediate-time regime. The results of the present 
paper will eventually constitute an essential ingredient in the construction of the atomic 
version of the self-consistent GLE theory. The most immediate task, however, is to provide 
a more formal understanding of the physical meaning and the fundamental nature of the 
kinetic friction represented by (°, and in general, of the Langevin equation for atomic liquids 
in Eq. (I2.2p . This task is addressed in the following two sections. 

III. KINETIC FRICTION ON AN ATOMIC TRACER PARTICLE. 

The fact that the Langevin equation in Eq. (12.21) adequately describes the results for 
the MSD of Newtonian systems strongly suggests a more profound equivalence between the 
dynamics of Newtonian and Brownian liquids, at least regarding tracer diffusion phenomena. 
In both cases the relaxation time t$ = [M/(°] of the velocity, due to the friction force 
— C°v(t), defines the crossover from ballistic (t <C t$) to diffusive (t ^> t$) motion. The 
fundamental difference lies in the physical origin of the friction force — (°v(t): in a Brownian 
liquid, this friction is assumed to be caused by an external material agent, namely, the 
supporting solvent, which also acts as a heat reservoir. In a Newtonian liquid, in contrast, 
the ('Doppler') friction force — C°v(t) is not caused by any external material agent but by 
the molecular collisions responsible to impose the equipartition of energy among the kinetic 
energy degrees of freedom of the system. 

In other words, the underlying physical origin of the kinetic friction effects is that for times 
t much longer than To, in which each particle has undergone many molecular collisions, the 
kinetic energy has indeed been distributed according to the most probable distribution. This 
then means that partial thermal equilibrium, consisting of this equipartition of the kinetic 
energy, is achieved within the time scale represented by r . Such partial thermalization must 
involve the transport of heat through molecular collisions. The instantaneous fluctuations of 
this local temperature, however, lead to the random force f°(t) that cause the instantaneous 



12 



fluctuations in the velocity of the tracer particle. The same molecular mechanisms are also 



responsible for the emergence of the systematic friction force — C°v(t) on the tracer particle 



whenever this particle has a non-zero mean velocity v(t). 

Let us now develop a more microscopic and precise discussion of these physical effects. 
For this, let us now consider (N + 1) particles in a volume V, the first of them (the tracer 
particle) having momentum p(t) = Mv(t) and the rest iV particles having momentum 
Pi(t) = mvi(t), with 1 < i < N. In the absence of external forces, the total momentum 
P total = p(t) + Y^Li P*(*) is conserved, 

rl N 

-[p(t) + ^ Pl (t)] = 0, (3.1) 

or 

-[p(t) + j tir,t)d 3 r] = 0, (3.2) 

with 

N 

j(r,*) = ^p^r-r^)), (3.3) 



«=i 



and with V' = V — Vt being the total confining volume V minus the excluded volume Vr 
around the center of the tracer particle. Thus, V is bounded by an outer surface £ of the 
confining walls, which remain fixed in space, and by the surface a around V, which follows 
the motion of this particle, and hence, is not fixed in space. Clearly, then, the momentum 
conservation equation can also be written as 

dp(t) f fdj(r,t) \ , 



m j v , v m ) dr - (3 ' 4) 

On the other hand, taking the time derivative of Eq. (13. 3p . one can write the continuity 
equation for the momentum density j(r, t) as 



Of 



V-nx(r,£) + V-IT/(r,£), (3.5) 



T7 TT 

where YIk ( r , t) and lit/ ( r 5 1) are the kinetic and configurational components of the stress 
tensor n (r, t), i.e., 

n(r,£)=n^(r,£)+fh7(r,£). (3.6) 
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The kinetic part TIk (r, t), which describes the change of momentum due to particles crossing 



the boundaries of V, is a second rank tensor with components 32] 

K K{r , t f^_±mm s(I _ im , ( 3. 7 ) 

L J *-^ m 

The configurational part of II (r,£), on the other hand, is defined by the condition 



N 



V- Ibv (r, t) = J2 MW(r ~ r;(*)), (3.8) 



i=i 



so that the continuity equation for j(r, £), Eq. (13. 5p . can actually be written most simply as 

Thus, substituting this equation in Eq. ([3.4]) . we have 

^ = -^(v-n^(r,t)) d 3 r-^ fj]p i (t)5(r-r i (i))J d 3 r. (3.10) 

The second term on the right side of Eq. ( I3.10J) is just — Ei=i Pi(t) = ~ Si=i ^i(t) = 
-Ef=iE^F,,(t) + F lT (t)] = -EiiF, T (t) = EtiF«(t), where F tJ (t) is the force 
exerted on particle % by particle j and Fxi(t) is the force exerted by particle i on the tracer 
particle, and where we have used the fact that Fy(t) = — F^(t). We may then write this 
term as £ii F„(t) = Eii V«ri(| r, - r T |) = /[V r «(| r - r T (t) |)] Eii *(r - r 4 (t)) rf 3 r. 
By shifting the origin of the coordinate system to the center of the tracer particle (including 
the variable of integration r), it is not difficult to see that Eq. (I3.10p can be rewritten as 

dY>{t) - j (V- U K (r,t)) d 3 r + ! [Vu{r)]n*{r,t) d\ (3.11) 

N 
i=l 

being the local particle density around the tracer particle described, however, from a refer- 
ence frame whose origin moves together with the center of this particle. 

Eq. ( 13. lip shows that there is a very simple and exact coupling between the force on 
the tracer particle and two collective variables of the surrounding fluid, namely, the kinetic 
component U.k ( r , t) of its stress tensor (whose trace is a measure of the instantaneous 

14 



dt 

with 



local kinetic energy per particle, i.e., of the local instantaneous temperature) and the local 
number density n*(r,t). Thus, the interatomic forces affect the motion of a tracer particle 
through two physically distinct channels. The first of them is a kinetic constraint imposed by 
momentum conservation, and the second is the ordinary configurational effect of interatomic 
forces. Each of these two variables contribute additively and linearly to the total force on 
the tracer particle. 

Taking the equilibrium average (indicated by an overbar) of Eq. (j3.1ip . we have that 



dp 
~dt 



- (V-rbr(r)J d 3 r + [Vu(r)]n*(r) d 3 r, (3.13) 

where the mean value of ELr- (r, t) can be obtained averaging Eq. (13. 7p . with the result 



11* (r) = -(k B T) I n*(r), (3.14) 



with I being the 3x3 cartesian unit tensor and with n*(r) given by 



n*(r) = ng(r), (3.15) 

where g(r) is the bulk radial distribution function of the system. From these results it is 
not difficult to realize that for symmetry reasons each of the two mean forces on the right 
side of Eq. (13.131) vanish independently, so that the tracer particle experiences a vanishing 
total mean force, and (dp/dt) = 0. 

We may now write the state variables IIk (r) and n*(r) as the sum of their equilibrium 
mean value plus the corresponding fluctuations, namely, as 



UK(r) = UK(r) + 5U K (r) (3.16) 

and 



n*(r) =n*(r) + 5n*(r). (3.17) 

This allows us to rewrite Eq. ( 13. lip as an exact relationship between these thermal fluctu- 
ations and the instantaneous momentum of the particle, namely 

dp(t) 



dt 



f (V-<*n*M)) d 3 r+ J [Vu{r))8n*{r,t) d 3 r. (3.18) 



This exact equation will now be taken as the starting point for a formal statistical mechanical 
derivation of the ordinary Langevin equation of a the tracer particle in our atomic liquid. 
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For this, let us first recall that the basis of the GLE formalism are the general mathe- 

matical conditions stated by the theorem of stationarity [12]. This theorem states that the 

equation describing the dynamics of the thermal fluctuations 5a,i(t) (= di(t) — af 7 ) of the 

instantaneous value of the macroscopic variables aiif) (i — 1, 2, ..., v) around its equilibrium 

value af 7 must have the structure of the most general linear stochastic equation with additive 

noise for the vector 5&(t) = [Sai(t), Sa^it), ..., 6a„(t)] , namely, 

t 
dSa(t) 



dt 



-cox^S&it) - I L{t- t'^Sa^dt' + f (£). (3.19) 



In this equation \ is ^ ne niatrix of static correlations, Xij — (^Oi(O)^a^(O)), to is an anti- 
Hermitian matrix (uij = —00%), and the matrix L(t) is determined by the fluctuation- 
dissipation relation L^it) = \fi(t)fj(0)), where fi(t) is the ith component of the vector of 
random forces f (£). Besides the selection rules imposed by these symmetry properties of the 
matrices y, oj, and L(t), other selection rules are imposed by other symmetry conditions. For 
example [12], if the variables ai{t) have a definite parity upon time reversal, a$(— t) = \idi(t) 
with Aj = 1 or -1, then u>ij = —XiXjUij and Lij(t) = \i\jLij(t). 

Let us now apply this mathematical infrastructure to the physical context involving the 
exact momentum conservation equation, Eq. ( 13. lip , and let us define the vector 5&(t), 
partitioned as 

5a(t) = [p(£), 5U(t), 5n*(t)f , (3.20) 

in terms of the sub- vectors p(t), SH(t), and 5n*(t), defined by their components 

p(t) = ( Px (t),p y (t), Pz (t)), (3.21) 

[(ffl(£)] Q/3 (r) = 8Uf(r, t) (with a,fi = x, y, z, and r e V) (3.22) 

and 

[5n*(t)}(r) = 5n*(r,t) (with r e V). (3.23) 

With this definition of the vector 5&(t) one can calculate the static correlation matrix 
X = (<5a(0)5a"l~(0)) using the microscopic definitions of 5n*(t) in Eqs. (I3.12p . (j3.15p . and 
(I3.17P and of SH^ (r, t) in Eqs. (13. 7p . (I3.14p . and (j3.16p . The result for \ can be written as 
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the following partitioned matrix 



Xpp o 

" Xnn Xnr 

whose sub-matrices have elements defined as [Xpp] a/3 — {p a P l3 )i [Xtttt] u ( r , r ') = 
(6ng(r)m%(i>)), [x.n] Q/3 (r,r') = (6nf(r)6n*(r>)), [Xn.]^(r,r') = (6n*(r)6ll%(r>)), and 
[Xnn]( r i r ') — (Sn*(r)Sn*(r')), given, respectively, by 

[ Xpp \ a P = (Mk B T)5 a p, (3.25) 

[x.«T^ u (r, r') = (k B T) 2 {[5 a p5, u + 5 a ,5^ + 5 au 5^] X s(r, r') + 5^5, uXd (r, r')}, (3.26) 
[ X _]^(r, r') = {xn,] aP (r, r') = -(k B T)5^ Xd (r, r'), (3.27) 



and 

[Xnn](r, r') = Xs(r, r') + Xd(r, r') 

where the self and the distinct parts of Xnn are defined as 



X s (r, r') = n*(r)S(r — r') 



(3.28) 



(3.29) 



and 



Xd(r, r ) = n*(r)n*(r) — n*(r) n*(r) (3.30) 

We then write up the generalized Langevin equation for our vector 8a(t) in the format 
of Eq. (I3.19p . For this, we first notice that all the variables, except p(£), are even func- 
tions under time-reversal. According to Onsager's reciprocity relations, and the general 
ant i- hermit icity of to and hermiticity of Lit) 12], we have that the only possibly non-zero 
sub matrices of u and L(t) are 



LO- 







pn 





- 


^p-K 


Upn 















and 



L(t) 



(3.31) 



L pp {t) 

L n7T [t) L nn [t) 
LUt) L nn {t) 



(3.32) 
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The determination of some of the non-zero elements of uj and L(t) is rather straightfor- 
ward. Thus the previous selection rules, along with the general format imposed by the GLE 
equation (13.191) . allows us to write the time-evolution equation for the sub- vector p(t) as 

— — = - [Upnix' 1 )™ + Upnix^mr] SIl(t) ~ [uJ p7T (X~ l )nn + UpniX^nn] Sfl* (t) 

- / dt'L pp (t - t'Xx-^ppPit') + f p (t) (3.33) 

Jo 

By comparing with the exact momentum conservation equation in Eq. (13.331) . we immedi- 
ately conclude that L pp (t) = f p (t) = 0, and that the remaining terms correspond, respec- 
tively, to the kinetic and configurational forces on the right side of this equation. In addition, 



for simplicity we approximate n*(r)n*(r) fa n*(r) n*(r) in Eqs. ( 13. 27ft and ( 13. 30 h . so as to 
neglect at this point the static cross-correlation Xnn and Xnn, 

[ X7rn ]^(r,r') = [Xn 7r ] a/3 (r,r')^0, (3.34) 

so that the previous equation is rewritten as 

— — — = -W plr (x~ 1 )iTir8'n-(t) ~ ^pn{X~ l )nnSn*{t). (3.35) 

Comparing this equation with Eq. H3.33|) one can determine the sub-matrices u pn and co pn 
and hence, also the sub- matrices oj np (= — wL.) and u np (= — u^J. 

In a similar manner, from the exact format imposed by the GLE, and using the previous 
selection rules (as well as the approximation Xnn ~ 0), one can also write the time-evolution 
equations for the other two variables, 5U(t) and Sn*(t), as 

^^ = -oo npXpp p(t) - J dt'L^(t-t')( X ~ l \JW) 

dt'L m {t - t'){x' l )nrM*{t') + U{t) (3.36) 



and 



dSn*(t) 



w«pX^P(*) - / dt'L n „(t - t')( X '%JW) 

dt'L nn (t - t')( X ~ l )nJn*(t') + f n (t) (3.37) 



'o 

Eqs. (I3.35p - ()3.37p provide a non-contracted description of the thermal fluctuations in an 
atomic liquid, which involves the tracer particle's momentum p(t) as one of the variables. 
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Contracting this description to the state subspace spanned by p(t) itself will finally lead to 
the complete generalized Langevin equation for a tracer particle in such atomic liquid. The 
result, however, is rather involved, but the essence may be best appreciated if we introduce 
an additional simplification, which consists of neglecting the dissipative coupling between the 
variables 5H(t) and Sn*(t), i.e., by setting L nn (t) = L nn (t) = 0. Under these circumstances, 
the solution of Eq. (13.361) can be written as 

5U(t) = G„(t)6Il(0)- [ dt'G^t-t'^pX^pit') 

Jo 

+ f dt'G^{t-t')U{t') 
Jo 

where the Green's function G n (t) is the solution of 

dG n {t) rt 



(3.38) 



dt 



[ dt'L^it-t'Kx-^Ait'), 

Jo 



(3.39) 



with initial condition G n (t = 0) = J, i.e., its Laplace transform (LT) G n (z) will be given, in 

terms of the LT of L nn (t), by 

-i 



GJz) 



zl + L n7T (z)(x l ) 



(3.40) 



Substituting the expression above for 5U(t) in Eq. (13.351) we arrive at the following 
Langevin equation 



dp(t) 
dt 



1 



M 



— dt' ( K (t - t') ■ p(t') + f K (t) - Upnix'XJn^t), 



where we have defined the time-dependent kinetic friction coefficient Cx(^) as 

•o- 

C (t) 

— = —U pn (x )mrG 7T (t)ljJ np X p p 



(3.41) 



(3.42) 



(3.43) 



According to the theorem of contractions [12j|, ix(t) and ( K (t) must satisfy the fluctuation- 
dissipation relationship (fj£ (t)f^(t')) = k B T(^f{t — if). 

For future reference, let us notice that Eq. (I3.40p allows us to write the FT of the kinetic 
friction function ( K (t) in Eq. (I3.42[) directly in terms of the FT of the memory function 
L irn [t) as 



and the kinetic random force fx(^) as 
ffir(t) = -Upnix' 1 ) 



G 7r (t)5U(0) + [ dt'G n (t-t')U(f) 
Jo 



M 



= -Upn(X ^tttt Zl + L^(Z){X l ) 



UivpX 



pp ' 



(3.44) 
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Clearly, the exact determination of the memory function L nn (t), and hence, of the kinetic 
friction coefficient ( K (£), is perhaps impossible, but some properties can be drawn from the 
expressions just derived, at least in certain limits and within well-defined approximations. 
For example, if ( K (t) relaxes to zero within a finite relaxation time, then for times t much 
longer than such relaxation time we can approximate ( K (t) by its Markov limit, 






( K (t)n25(t)C, (3.45) 

where 



C°= / dt( K (t). (3.46) 

Jo 

<-> 
In addition, due to the radial symmetry of the interparticle interactions, (° must be isotropic 

(i.e., diagonal), so that 

f°= C° I, (3-47) 



with (° given by 



C° = -M 



UJ 



L-l(z = Q)^ pXp i ' • (3-48) 



Thus, we conclude that in the Markov limit Eq. (13.411) can be written as 

M^ = -C°v(t) + f°(t) - uj pn {x' l )nJn\t). (3.49) 

In the following section we discuss the additional contraction process that leads to the 
elimination of the variable Sn*(t) from the description consisting of this equation and of 
Eq. (I3.37P above. Before that, however, let us mention that Eq. ( 13.441) . and in particular 
its Markov limit in Eq. (13. 48 p . involves the LT of the memory function L~£(z) as the only 
unknown quantity, which the GLE formalism is unable to determine. Although one could 
introduce additional approximations to determine this memory function, this is not the main 
objective of the present paper; instead, the derivation above was only meant to provide a 
more formal explanation of the origin of the kinetic friction and random forces, introduced 
and discussed in more efficient and intuitive terms in the previous section. After all, such 
arguments did provide a simple and accurate zeroth-order approximate determination of the 
kinetic friction coefficient (°, namely, the use of the kinetic-theory value of the self-diffusion 
coefficient, Eq. ( 12. 4(1 . in Einstein's relation, Eq. (12. 3ft . 
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IV. CONFIGURATIONAL FRICTION ON AN ATOMIC TRACER PARTICLE. 

One important contribution of the previous section was to make a point that in an atomic 
liquid the force on a tracer particle couples linearly with the kinetic component Hk (r, t) 
of the stress tensor and with the local number density n*(r,t) of the surrounding fluid, 
as indicated by Eq. (13. lip . As a consequence, the interatomic forces affect the motion of 
a tracer particle through two physically distinct channels, namely, the kinetic constraint 
imposed by momentum conservation and the ordinary configurational effect of interatomic 
forces. Another important conclusion was to notice that the former is the origin of the 
kinetic friction, finally formatted in Eq. (I3.49P as a dissipative friction term — C°v(t) plus 
the corresponding random force f°(t). Let us now discuss the effects of the coupling with 
6n*(r,t). 

For this, let us resume the formal process of contraction of the description initiated in 
the previous section. We thus recall that after projecting out the variable 6 Uk (r, t), the 
time-evolution equations for the remaining fluctuating variables are eqs. (13.411) and (13.371) . 
For clarity, we rewrite here these equations as 

^ = -^P(«) + f °(*) - WpnOr 1 )™^), (4.1) 

and 

^P- = -u np x£v{t) - J dt'L nn (t - t'^x-XJn^t') + f n (t), (4.2) 

where the kinetic friction term of eq. ( I3.4ip has been written in its markov limit (as in Eq. 
(I3.49p ). and where the term of eq. (13.371) involving L n7T (t) has been neglected, as discussed 
immediately above Eq. (13.381) . We now formally project out the variable 8n*(t) by solving 
Eq. (14. 2 p for this variable, and substituting the resulting solution in the third term of the 
right side of Eq. ( 14. ip . 

This contraction process results in the following generalized Langevin equation for the 
velocity v(t) (= p(t)/M) of a tracer particle in the atomic liquid, 

M^- = -C°v(t) + f°(t) - |* dt> AC(t -f) ■ v(0 + F(t), (4.3) 

where the configurational time-dependent friction function A((t) is given by 

M = -U pn ( X - l )nnG(t)u npX -^ (4.4) 
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and the new (configurational) random force F(t) as 



" [v ) — ^pn\A )nn 



G(t)6n*(0)+ I dt'G(t-t')f n (t') 
o 



(4.5) 



In these equations, the Green's function G(t) is the solution of 

dG{t) rt 



dt 



dt'LrUt-tKx-XnGit') (4.6) 



with initial condition G(0) = I. According to the contraction theorem [12|, F(t) and A ( (t) 
must satisfy the fluctuation-dissipation relationship (F a (t)F /3 (t')) = k B TAC al3 {t — £'). 

Although the previous statements are physically accurate and well-defined, it is also useful 
to rephrase this abstract derivation in more concrete and intuitive terms. For this we rewrite 
Eq. (jUD as 

M rfv(t) = _ c o v(t) + f o (t) + J d 3 r [Vw(r)]5n*(r, t), (4.7) 

to recover the original notation in Eq. (J3.18P for the configurational force term. The com- 
parison with Eq. (14.11) then implies that [— w pn (x~ 1 )„„] a (r) = [V a w(r)] and, convoluting 
this equation with [% nn ](r, r') = (5n*(r)5n*(r')) determines that u pn is given by 

K„] a (r) = - I dV[V' a M (r')]Xnn(r\ r) 

= k B T[V a n eq (r)l (4.8) 

where the second equality is a direct consequence of the exact equilibrium condition referred 
to as the Wertheim-Lovett's relation, namely [33l |. 



[W(r)] = -pJd 3 r' X nn(r',r)[V a u(r')}- 



(4.9) 



In the previous equations, the equilibrium mean value n(r) has been denoted by n eq (r). 

Now, since u np = —[u pn W we find that [u; n p] a (r) = — k B T[V a n eq (r)}. Using this result, 
along with the value [Xpp] a/3 = (Mk B T)5 a i3 (Eq. (13.251) ). we can write Eq. ( 14. 2 p more 
concretely as 

d6n *j r >^ = [\/n eq (r)} ■ v(t) - [ dil [ d 3 r'D*(r, r'; t - t')5n*(r', t') + /(r, t), (4.10) 
dt J J 

with -D*(r, r';t) being the elements of the "matrix" D*(t) = [L nn (t)(x~ 1 ) n n\- The first term 
on the right side of this equation is a linearized streaming term and f{r,t) is a fluctuating 
term, related to D*(r, r'; t) by (f(r, t)f(v', t')} = J d 3 r"D*(v, r"; t - t')xnn(?" , r')- 
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Formally solving Eq. f |4.10|) and substituting the solution for 5n*(r, t) in Eq. (14.71) . leads 
again to the generalized Langevin equation in Eq. (14. 3 j) with the time-dependent friction 

■H- 

tensor A£(£) of Eq. (14.41) given more concretely by 

A((t) = - j d 3 r f d 3 r'[Vu(r)]G*(r,r';t)[Vn eq (r% (4.11) 

where G*(r, r';£) is the propagator, or Green's function, of Eq. (14.101) . i.e., it solves the 
equation 

d^'for';*) =- f dt' f d 3 r"D*(r,r";t -t')G*(r",r';t'), (4.12) 

with initial value G*(r, r'; t — 0) = <5(r — r'). Notice that, since the initial value Sn*(r, t = 0) 
is statistically independent of v(t) and f{r,t), the density-density time-correlation function 
X*(r, r'; t) = (5n*(r, t)Sn*(r', 0)), which is the van Hove function of the particles surrounding 
the tracer particle, and observed from the tracer particle's reference frame, is also a solution 
of the same equation with initial value x*( r ) r 'i t = 0) = Xnn{r, **')■ 

To simplify the notation, let us re- write Eq. <|4TTT|> as A C (t) = -[Vu f ] ■ G*(t) ■ [Vn e «], 
where the convolution J d 3 r"A(r, r")B(r", r') between two arbitrary functions A and B is 
written as the inner product A ■ B, and similarly with (column) "vectors" such as u and 
n eq . In this notation, the dagger means transpose. With this notation, Wertheim-Lovett's 
relation reads [Vn e<? ] = —ftXnn • [Vw] . With this relation, and the definition of the the inverse 
matrix Xnn by the equation XnnXnn = I, with I being the unit matrix (J(r, r') = <5(r — r') 
= Dirac's delta function), one can write Eq. (14. lip in a variety of different but equivalent 
manners. In particular, we will employ the following: 

A C (t) = k B T[Vn^\ ■ x~nl ■ X*(t) ■ Xnl ■ \^rf\ (4.13) 

where we have used the fact that the van Hove function x*{t) can be written as x*{t) = 

G*(t)- X nn. 

Let us now notice that for spherical particles A ( (t) is isotropic and diagonal, 

AC(t)=/AC(t), (4.14) 

so that we only have to calculate the scalar time-dependent friction function A((t). The 
exact expressions for A£(£) can then be given a more concrete and tractable appearance 
if some approximations are introduced, related to the general properties of the functions 
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X*(r, r';t) and Xnn(*, r'). The latter is just the two-particle distribution function of the 
colloidal particles surrounding the tracer particle, but subjected to the "external" field u(r) 
exerted by this tracer particle. Thus, it is effectively a three-particle correlation function. 
Only if one ignores the effects of such "external" field, one can write x nn (r, r') = Xnn{\ r — 
r'|) = n5(r — r') + n 2 [g(|r — r'|) — 1]. Similarly, we may also approximate y*( r > r';£) by 
X*(|r — r'|; t). This is referred to as the "homogeneous fluid approximation" [ll|, which then 



allows us to write 



and 



with 



and 



Xnn(r,r';t) = (1/2tt) 3 f d 3 kexp[ik ■ r]nS(k) 

X*(r,r';£) = (1/2tt) 3 / d 3 kexp[ik ■ r]nF*(k, t), 

1 - 
F*(k,t) = -(Y,ew\ik-[v t (t)- rj (0)]]). 



i, j 



(4.15) 
(4.16) 
(4.17) 

(4.18) 



S(k) = F*(k,t = 0). 

Using these expressions in Eq. (14.131) . along with the fact that Vn e? (r) = nVg(r) = riWh{r) 
(so that its FT is ihnh(k) = ik.[S(k) — 1]), we have that Eq. (14 . 1 3 j) becomes 



AC(t) 



k B T 



3(2tt) n 



dk 



k[S{k) - 1] 



F*(k,t). 



(4.19) 



The function F*(k,t) in this equation is just the intermediate scattering function, but 
the asterisk indicates that the position vectors rj(£) and i"j(0) have the origin in the center 
of the tracer particle. Denoting by xy(£) the position of the tracer particle referred to a 
laboratory-fixed reference frame, we may re-write 

"l N 
-JjexpCzk-teW-x^O)]) 



F*(k,t) = { 



'•.y 



[exp(ik.[x T (*)-x T (0)])]), 



(4.20) 



where rj(t) is the position of the ith particle in the fixed reference frame. Approximating 
the average of the product in this expression by the product of the averages, leads to 



F*(k,t)^F(k,t)F s (k,t), 



(4.21) 



where Fs(k,t) = (expjik • \x r(t) — xr(0)]}) is the self ISF. This is referred to as the 
decoupling approximation [ll| . Thus, from the exact result in Eq. ( I4.13P above, plus the 
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introduction of the two approximations just described, we finally arrive at the following 
general but approximate expression for the time-dependent friction function A£(t), 

'k[S(k)-l]' " 



A((t) = kB \ I dk 



F(k,t)F s (k,t). (4.22) 



S(k) 

This expression is reminiscent of the corresponding mode coupling theory (MCT) result 
9|. Its derivation above, however, follows a completely different conceptual route. In the 
following section we discuss important implications of our results above. 

V. LONG-TIME DYNAMIC EQUIVALENCE. 

The main result of the previous sections is, of course, the generalized Langevin equation 
(Eq. (14.31) ) describing the ballistic to diffusive crossover of the Brownian motion of individual 
tracer particles in an atomic liquid. Taking into account the isotropy of the configurational 
time-dependent friction function (Eq. (14.14)) ). this stochastic equation reads 

M^- = -C°v(t) + f°(t) - f dt'A({t - i')v(t') + F(t). (5.1) 

dt Jo 

The configurational effects of the interparticle interactions is embodied in the time-dependent 
friction function A((t), which Eq. (14.22)) writes in terms of the ISFs F(k,t) and Fs(k,t). 
Thus, the full analysis of this stochastic equation requires in principle the previous determi- 
nation of these more complex dynamic properties. Some important implications, however, 
can be drawn without a detailed knowledge of A£(£). 

The most remarkable conclusion is that the Brownian motion of individual tracer particles 
in atomic and colloidal liquids is described by the same mathematical model, namely, the 
generalized Langevin equation derived here for atomic systems (Eq. (15. ip with Eq. (14.22)) ). 
and the GLE derived in Ref. [llj for colloidal fluids (Eq. (Oj) with Eq. (TOD ). According 
to this formal mathematical similarity, the properties that describe the tracer's random 
motion in atomic and in colloidal liquids, such as the mean squared displacement W(t), 
should collapse onto each other when expressed in dimensionless units that absorb the mass 
M and the short-time friction coefficient (Q s > or £°). 

To see this, let us first notice that from Eq. ( 15. ip one can write the velocity autocorrelation 
function V(t) in terms of A((t), in Laplace space, as 

k B T 

v{z) = TTsn-mi ' < 5 - 2 > 

c ' M ~ M 
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with the friction coefficient (s representing either the kinetic friction coefficient £° in atomic 
fluids or the solvent friction coefficient (^ in colloidal liquids, 



r (° (for atomic systems) 
£( s ) (f or Brownian systems) 



Cs = { ,., _ (5.3) 



Using this result in the exact relationship in Eq. (12.1 Op one can derive the following integro- 
differential equation for W(t), 



at j 



AC(t - g) 

(s 



W{t')dt', (5.4) 



where 

rs = M/Cs (5.5) 

is the crossover timescale from ballistic to diffusive motion (and which equals the mean 
free time To only in atomic liquids), and where the short-time self-diffusion coefficient D$ is 
defined by Einstein's relation, 

D s = k B T/( s . (5.6) 

Thus, Ds = D° is given by the kinetic theoretical result in Eq. ( 12. 4p only for atomic liquids. 

For both, atomic and Brownian tracers, in the absence of interactions A((t) vanishes and 
Eq. ( 15. 4p becomes Eq. ( 12. lip , discussed in Sec. |TT]for atomic liquids, and whose solution is 
given by Eq. (12.12H . In the presence of interactions, however, A((t) ^ 0, but at very short 
times (t ^C Ts) the solution of Eq. (15. 4p is still identical to that of a freely-flying particle, 
i.e., the short-time asymptotic expression for the MSD is also given by W(t) w ^o^ 2 ? as 
illustrated by the molecular dynamics simulations in Fig. [TJ 

In the opposite regime, t 3> Ts, the interparticle interactions change the long-time asymp- 
totic limit of W(t) from its free-diffusion value W(t) ~ D°t to the new value W(t) ~ Drf, 
which defines the long-time self-diffusion coefficient Dl. In this regime, the convolution 
j A((t — t')W(t')dt' on the right side Eq. (15. 4p can be approximated by its Markov limit 
[r o °° A((t')dt']W(t) = A(W(t), so that Eq. (E2D reads 

dW(t) W{t) 
rs^ r + ^ = D s t, (5.7) 



where 



D ~ Ds~ 1 + AC fts' ( 
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with the constant A( defined as A( = f Q °° A((t)dt. This equation is precisely Eq. ( 12. 18[) 
of Sec. HIl whose analytic solution in Eq. (I2.17P was shown in Fig. 1 to provide a simple 
interpolation between the short- and long-time limits of the molecular dynamics simulation 
data for the MSD of the HS liquid throughout its thermodynamically stable liquid regime, 
O<0<O.5. 

To continue the discussion of the dynamic equivalence between atomic and Brownian 
liquids, let us define in both cases a length Is in terms of t$ and D$ as 

l 2 s = D s r s . (5.9) 

For atomic liquids Is is identical to the mean free path / , but not for colloidal systems, for 
which Is is only given by this equation (together with Eqs. ( 15. 5 p and (I5.6P ). We may now 
use Is and Ts as the units of length and time, respectively, and rewrite Eq. (15 .4p in terms 
of the scaled time t* = t/r s and the scaled MSD w(t*) = W(t)/l s as 

dw ^ + w{t*) = t* - ( — ) [ A(*(t*-t*')w(t*')dt*', (5.10) 



dt* \Tj J Jo 

where the configurational timescale r/, the time it takes a particle to diffuse a mean distance 
d = n -1 / 3 with a diffusion coefficient Ds, is given by 

T! = d 2 /D s , (5.11) 

and in which we have defined the dimensionless function A£*(£) as 



AC(t) = 



r/AC(t) 



dk 



k[S(k) - 1] 
S(k) 



F(k,t)F s (k,t). (5.12) 



" 3(2vr) 3 n 5 /3 

Let us notice that the purpose of using lg and t$ as the units of length and time is to 
focuss on the crossover time-regime from ballistic to diffusive motion. If, instead, we were 
interested in focussing on the crossover from free- diffusion to structural relaxation, the best 
would be to use the mean inter-particle distance d = n -1 / 3 as the unit of length and 17 as 
the time unit, and to rewrite Eq. (j5.4p in terms of the scaled time t* = t/ri and the scaled 
MSD w(t*) = W(t)/d 2 , to read 

(—] < ^$^ + wit*) =t*- f A(*it* - f')w(t*')dt*' . (5.13) 

\TiJ dt* J 

For dense atomic liquids (e.g., hard spheres at <f> < 0.5) either choice is perfectly adequate to 
observe within the same time window both, the crossover from ballistic to diffusive motion 
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and the crossover from free to correlated motion, as illustrated by the simulation results for 
W(t) in Fig. 1. The reason is, of course, that in this case ts and 77 do not greatly differ 
from each other, and hence, there is no important time-scale separation. 

The situation is, however, dramatically different in the corresponding colloidal liquid 
since, as discussed in Sect. [Til the ratio ts/tj may be as small as t s /t i ps 1CT 8 . Thus, these 
two crossover timescales are separated by about 8 decades, and cannot be analyzed in the 
same time window. In fact, this implies that if we focus on the crossover from ballistic to 
diffusive motion, as in Eq. (j5.10p . then the term involving the configurational friction will 
be completely negligible. In contrast, if we focus on the crossover from free to correlated 
motion, as in Eq. ( I5.13p . then it is the inertial term involving the time derivative of the MSD 
what can be neglected. This is referred to as the overdamped limit, which also amounts to 
ignoring the inertial term M[dv(t)/dt] on the left side of the GLE in Eq. (j5.ip . Thus, in 
this limit Eq. (J5.13P reads 

w (t*) =t*- AC(t* - t*')w(t*)dt* . (5.14) 

Jo 

This only changes the true short-time limit w(t*) ~ (d/ls) 2 t* 2 /2 of the solution of Eq. ( 15. 131) 

to w(t*) ~ t*, but leaves unaltered the long-time limit w(t*) ~ D*t*. This equation thus 

describes the diffusive motion of colloidal tracer particles. 

It should also be clear, however, that even though in atomic liquids there is not an appre- 
ciable timescale separation, Eq. (j5.14p also describes the long-time motion of atomic tracer 
particles. Hence, except for the referred short-time differences, the MSD of an atomic and 
a colloidal liquid with the same interactions and the same S(k) should be indistinguishable 
when plotted in terms of these dimensionless units, provided that the intermediate scatter- 
ing functions F(k,t) and Fs(k,t), which enter in the previous expression for A£*(£) above, 
also share a similar long-time scaling property. Thus, we may embark on a study of these 
dynamic properties to see if they indeed exhibit the desired scalings, or else, we can check 
directly if the MSD itself exhibits the expected universality of atomic and colloidal liquids. 

Here we have adopted the second approach, and have compared the molecular and the 
Brownian dynamics simulation data of W(t) for the HS system at the volume fractions 
cf) = 0.1 and 0.5. For this comparison we use the same molecular dynamics data as in Fig. (U 
now plotted in terms of the dimensionless MSD [H / (t)/6c? 2 ] as a function of the dimensionless 
time tjri = [D°t/d 2 ] , and the Brownian dynamics data generated for this comparison using 
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FIG. 2: Dimensionless mean squared displacement [IF(i)/<i 2 ] as a function of the dimensionless 
time [Dst/d 2 ] for the volume fractions <fi = 0.1 and 0.5. The squares correspond to Brownian 



Dynamics simulations g, 



34l ] generated as described in Ref. |31|. and the circles to the molecular 



dynamics simulations of Fig. 1. 



the methodology explained in Ref. 3l|. As we can see from this comparison, for each volume 
fraction the molecular dynamics and the Brownian dynamics data agree at long times, within 
a high degree of numerical precision. The short-time difference between the molecular and 
Brownian dynamics simulation data originates, of course, from the fact that the latter are 
based on the conventional Ermak and McCammon's Brownian dynamics algorithm 8|, |34| . 
in which the "overdamped" limit is previously taken in the microscopic equations of motion. 
Thus, this agreement is also an indirect indication that the intermediate scattering functions 
F(k,t) and Fs(k,t) must also share similar scaling properties. The analysis of this issue, 
however, is addressed separately [35 ] . 



VI. DISCUSSION AND SUMMARY. 



In summary, in this paper we have explained a simple argument that exposes a dynamic 
equivalence between the long-time dynamic properties of atomic and colloidal liquids. Such 
simple arguments were complemented by a more formal fundamental derivation of the gen- 
eralized Langevin equation for a tracer particle in an atomic liquid, which is the atomic 
counterpart of the GLE for tracer diffusion derived in Ref. [llj for colloidal liquids in the 
absence of hydrodynamic interactions. The dynamic equivalence suggested by the fact that 
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the GLE for tracer diffusion in both cases has the same mathematical structure, need in 
reality that other dynamic properties, such as the intermediate scattering functions F(k, t) 
and Fs(k,t), also share a similar long-time scaling property. For the time being, here we 
have tested the predicted dynamic equivalence at the level of the mean squared displace- 
ment W(t) in the context of a specific model system, namely, the hard sphere liquid, in its 
dynamic version corresponding to molecular and Brownian dynamics. 

Let us state, however, that the present work does not settle the question of the generality 
of this dynamic equivalence. Instead, it only contributes to stimulate the corresponding 
discussion. For example, it is important to discuss the manifestation of this dynamic equiv- 
alence on properties other than the MSD. As indicated above, verifying that similar scalings 
are exhibited by the intermediate scattering functions F(k,t) and Fg(k,t) is an issue that 
must still be addressed in detail. In fact, our group has already approached this issue within 
the GLE formalism 35j, and the results turn out to be completely consistent with those 
of the present paper. The atomic-to-Brownian long-time dynamic equivalence thus seems 
to be a very robust prediction. The most relevant implications of this dynamic equivalence 
have been corroborated by the systematic comparisons between molecular and Brownian 
dynamics simulations of the sort illustrated in this paper. A summary of this analysis has 
been advanced in a recent brief communication [36( . Another important issue refers to the 
actual universality of the atomic-colloidal dynamic equivalence discussed in this paper, since 
the only validation of these predictions involved a specific model system, namely, the hard 
sphere fluid. 

In this direction, let us mention that the present colloidal-atomic dynamic correspon- 
dence is not restricted to the hard-sphere fluid, but it actually extends over to systems with 
soft repulsive interactions. This is a direct result of combining the present colloidal-atomic 
correspondence for the hard sphere system, with another important scaling rule, which de- 
rives from the principle of dynamic equivalence between soft-sphere and hard-sphere liquids 
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37| . The extension of this scaling to atomic systems is immediate once the collision 



diameter a entering in the expression for D°(n,T) in Eq. ( 12 Ah is given a proper definition 



381 ] for the soft-sphere potential u(r) considered. 

Still another issue refers to the possible limitations of this long-time dynamic equiva- 
lence, imposed by the fact that the present derivation apparently assumed colloidal systems 
in the absence of hydrodynamic interactions. In reality, however, the validity of this dy- 
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namic equivalence should extend over to systems with hydrodynamic interactions, provided 
that the corresponding effects enter only through the value of the short-time self-diffusion 
coefficient D s , as suggested in 39]. Other interesting directions along which to question the 
applicability and universality of this dynamic equivalence refers to the realm of liquid mix- 
tures and to the effects of attractive interactions. The answer to these questions, however, 
will only come from the comparison between the dynamic properties of atomic and colloidal 
liquids, similar to that presented here in Fig. 2, or by expanding the theoretical analysis 
that led us to the present proposal. 
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